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Abstract. We use the Gross-Pitaevskii equation to determine the spatial structure of the condensate density 
of interacting bosons whose energy dispersion e& has two degenerate minima at finite wave-vectors ±q. 
We show that in general the Fourier transform of the condensate density has finite amplitudes for all 
integer multiples of q. If the interaction is such that many Fourier components contribute, the Bose 
condensate is localized at the sites of a one-dimensional lattice with spacing 2ir/\q\; in this case Bose- 
Einstein condensation resembles the transition from a liquid to a crystalline solid. We use our results to 
investigate the spatial structure of the Bose condensate formed by magnons in thin films of ferromagnets 
with dipole-dipole interactions. 

PACS. 03.75.Hh Static properties of condensates; thermodynamical, statistical, and structural properties 
75.10.JmQuantized spin models - 75.30.Ds Spin waves 



1 Introduction 



This work is motivated by the recent discovery 
of a new coherence phenomenon of magnons in thin stripes 
made of the magnetic insulator yttrium- iron garnet (YIG). 
The energy dispersion e^ of the lowest magnon mode in 
this system has a rather unusual momentum-dependence 
which is crucial to understand the experiments: for a cer- 
tain range of orientations of an external magnetic field, 
efc exhibits two degenerate minima at finite wave-vectors 
+q and —q. The value of q is determined by a subtle 
interplay between exchange interactions, dipole-dipole in- 
teractions, and finite-size effects [MI]- The experimentally 
observed strong enhancement of the occupation of the 
magnon modes with wave- vectors ±q has been interpreted 
as Bose-Einstein condensation (BEC) of magnons. 

Another class of boson systems where the energy dis- 
persion has degenerate minima at finite wave-vectors ±qr 
are magnon gases in quantum helimagnets in a magnetic 
field 8J. Apart from the experimental relevance, the gen- 
eral problem of BEC in systems where the energy disper- 
sion has a minimum for non-zero wave-vectors is inter- 
esting on its own. In this context we mention the work 
by Yukalov [5] , who has investigated BEC in an interact- 
ing Bose system whose energy dispersion is minimal on a 
sphere in momentum space. He found that in this case the 
condensed state neither exhibits off-diagonal long-range 
nor is it superfluid. Moreover, Yukalov also pointed out an 
interesting analogy between BEC at finite momentum and 
the liquid-crystal phase transition, which can also be un- 
derstood in terms of a Ginzburg-Landau functional whose 



Gaussian term exhibits minima on a surface in momentum 
space PUTT]. 

The fact that BEC of quasi-particles is not necessar- 
ily accompanied by superfluidity has been emphasized by 
Kohn and Sherrington [T3] , who classified bosons into two 
different types: the first type consists of bound complexes 
of an even number of fermions; in the case of condensation 
of these bosons, superfluidity and off diagonal long-range 
order occurs. The second type of bosons consists of quasi- 
particles such as excitons and magnons; when the second 
type condenses, there is no superfluidity, but a change 
of spatial or magnetic order p~2].rT3] . Obviously, BEC of 
magnons in YIG is of the second type. From the point of 
view of critical phenomena it is not surprising that BEC 
at finite momentum is rather different from BEC at zero 
momentum. In fact, phase transitions which are character- 
ized by an order parameter which condenses on a surface 
in momentum space belong to their own universality class, 
the so-called Brazovskii universality class [THI15] . 

In this work we shall examine the general problem of 
BEC in a Bose gas whose energy dispersion has degenerate 
minima at two finite wave-vectors ±q. We show that in 
this case the time-independent Gross-Pitaevskii equation 
implies that the Fourier transform cf>k of the condensate 



wave-function has finite amplitudes 



?±q,(P±2q, ■ 



for 



integer multiples of the fundamental wave-vector q. Pre- 
viously a theoretical analysis of magnon-BEC in YIG has 
been performed by Tupitsyn, Stamp, and Burin [16] . How- 
ever, the effect of the spatial structure of the condensate 
wave-function as implied by the Gross-Pitaevskii equation 
was not considered by these authors. 
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2 BEC at finite momentum 

In this section we shall study BEC in a general class of 
interacting boson models on a lattice whose Hamiltonian 
is of the form 

H = H 2 + H 3 + H^ (1) 

The quadratic part H 2 of the Hamiltonian is given by 



H 2 = £ 



t , 7fc 



a W- k 



% 



fl-fcOfc 



(2) 



,t 



where a\~ and a k are the usual canonical annihilation and 
creation operators, the energy dispersion e^ is assumed 
to exhibit two degenerate minima at finite wave-vectors 
±q, and the terms proportional to the complex parame- 
ter 7fc explicitly break the U(l) symmetry associated with 
particle number conservation. In YIG these terms are re- 
lated to an external pumping filed, as explained in the 
appendix. In the absence of U(l) symmetry, the Hamilto- 
nian can also contain contributions involving three powers 
of boson operators, which in general are of the form 



H: 



1 



'< — Z^ fik 1 + k 2 +k a fi 

k 1 k 2 k 3 



1 



1 



j^aaa , j^aaa 

2 1 l;23 a -l a 2«3 + ^ 12;3 a -l a -2 a 3 

j^aaa , i p555 J t _t 

— 1 123 a 1 a 2 a 3 + —1 123 a_ 1 a_ 2 a_ 3 



(3) 



where N is the total number of sites of the underlying 
lattice. For simplicity we write a\ = a kl and abbreviate 
the interaction vertices by r k aa k k = r^fi etc-- Finally, 
the part H4 to the Hamiltonian involving four powers of 
the boson operators is in the absence of [/(l)-symmetry 
given by 



fci ...k4 



1 



(2!) 



i^aaaa t „t „ „ 
2 1 12;34 a -l a -2 a 3 a 4 



1 

3!' 



t J J 



7-iaaaa , j^aaaa 

■* I;234 a -l a 2fl3a4 + ^yi I23;4 a -l a -2 a -3 a 4 



ipaaaa „ „ „ i j^aaaa \ t t „t 

7[ J 1234 aia 2 a 3 a4 + —1 1234 a_]a_2a_3a_ 4 



(4) 



In the appendix we shall show how to obtain a boson 
Hamiltonian of the above form from an effective spin 
Hamiltonian describing the lowest magnon band of YIG 
in the so-called parallel pumping geometry 7.1MT?]. 

The spatial dependence of the Bose condensate is de- 
termined by the Gross-Pitaevskii equation [IHT , which can 
be obtained from the cxtremum of the corresponding Eu- 
clidean action. In order to write the various interaction 
processes in a compact notation, we introduce a two- 
component complex field <P k (r) , where r is the imaginary 
time and a — a, a labels the two components according to 
the prescription 



«fc -► n(r) 



J 



*i*(r). 



(5) 



The quadratic part H 2 of the Hamiltonian corresponds 
then to the Gaussian action 



S 2 [$\ = \f dr^i^-k^-k) 
1 JO u 



T + e fe - /i 

7fc 



Ik 



d T + efe - n 






(0) 



where /? is the inverse temperature and fi is the chemical 
potential. The Euclidean action corresponding to the in- 
teraction parts H 3 and H4 can be written in the following 
symmetrized form, 



L+fe 2 +fe 3 ,0 



" v kik 2 k 3 (J\<j 2 <j 3 

xir 3 (fc 1 a 1 ,fc 2f T2,fc3^Ki^ 2 2 ^3, 
rP 1 

s 4 [<p] = / dr— J2 J! 4 1+ ...+fc 4 .o 



h±...h^ <J\ ...<T4 



x T r A {k l o l ,...,k A ^W k \---$Z 



4! 



(7) 



(8) 



where the flavor indices <7i — a, a keep track of the two 
different field types, and the vertices r 3 {kiu\, k 2 a 2 , k 3 <j 3 ) 
and rjiffeio"!, . . . , ki(Xi) are completely symmetric under 
the permutation of all indices. The combinatorial factors 
in these expressions are chosen [TS] such that for a given 
ordering of the indices the completely symmetrized ver- 
tices can be identified with the partially symmetrized ver- 
tices appearing in Eqs. (|3l [4J) , for example 



r 3 (feio, k 2 a, k 3 a) = T{ 



;23- 



(9) 



In the presence of a Bose condensate some of the expecta- 
tion values (/)% — (4> k ) are finite and proportional to y/N, 
In equilibrium the order parameter fields 4> k are indepen- 
dent of the imaginary time. It is then useful to shift the 
integration variables <P in the Euclidean functional inte- 
gral according to ^%(t) = <t> k + 5<P k {T) an d expand the 
Euclidean action S[@\ = S 2 [<P] + S 3 {<P] + S^} in powers 
of the fluctuations, 



s[* + w = sw + f\ZjJifc 



sn(r) 



( 10 ) 

The physical order parameter field is determined by de- 
manding that the first variation of the action vanishes, 
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which yields the Gross-Pitaevskii equation 
6S[0\ 







&=<t> 






= (e fe -M)C fe +7feC fe 



k±k 2 "\oi 



Tf Z^ Z.^ 3k+k 1 +k 2 +k 3 ,o 



kik 2 k 3 o- 1 a 2 a 3 
1 



x J\ r ^ ka ' kl<71 ' fe2<J2 > ksa ^k\ < i>kl ( t ) 'k 



(ii) 



where we have defined 7^ = 7^ and 7^ = 7*,. 

To begin with, let us assume that the system condenses 
in a state where only the fc = mode is macroscopically 
occupied. Such a state tends to be favored if the dispersion 
has a minimum at k = 0. In this case 



3fc.O 



AW, 



(12) 



where the complex parameter ipQ = (V>o)* is expected 
to be of the order of unity. Assuming for simplicity that 
7o = 7o is real, we then obtain from our general Gross- 
Pitaevskii equation (11) 



= roV4 + 7o% +\J2 r3 ( 0<7 ' 0cri ' O^WW 

(71(72 

+ T\ £ A(0<r,0tr ls 0o-2,0o-3W^^, 



where 



^o = eo - M- 



(13) 



(14) 



We adopt here the standard notation in the field of criti- 
cal phenomena |20j where a negative value of tq implies a 
finite expectation value of the order parameter field. As- 
suming further that the three-legged vertices vanish and 
that only the particle number conserving four-legged ver- 
tex -Tot^oo = u 4 is finite and positive, we find that the 
Gross-Pitaevskii equation ( 13 1 has two degenerate solu- 
tions 



r = ±u 



'2(70 -r ) 



"4 



(15) 



Note that even for positive r = e — /i the condensate is 
stable if the energy scale 70 associated with explicit sym- 
metry breaking is sufficiently large. Of course, for 70 7^ 
there is no spontaneous symmetry breaking so that there 
is no gapless Goldstone mode in the condensed state. Due 
to the symmetry breaking terms a k a_ k and a_fcOfc in the 
quadratic part Hi of our Hamiltonian the effective poten- 
tial 

U eB [fa C] = N-'S^l -> \^V4,o^] (16) 

has two degenerate minima at purely imaginary values 
of the field as shown in Fig. [I] Cubic terms (which are 
neglected in these plots) distort the effective potential and 



r <0 




Retp 



r >0 




Im^ 



Reip 



Fig. 1. (Color online) Effective potential for zero momen- 
tum BEC in the presence of explicit symmetry breaking, see 
Eq. ( |16[ ). The quadratic part of the Hamiltonian is then given 
by Eq. |[2]|, where for simplicity we assume that 70 is real 
and positive. For the graphs the cubic vertices have been ne- 
glected and only the particle-number conserving component 
of the four-point vertices Foo^oo — u 4, has been retained. The 
graphs are for 70 — ro > where the effective potential has 
two degenerate minima on the imaginary axis. Upper graph: 
for ro < the center of the effective potential is a local maxi- 
mum so that its shape resembles Napoleon's hat. Lower graph: 
for ro > the local maximum in (a) transforms into a saddle 
point. 



break the degeneracy of the two minima. The problem of 
BEC at zero momentum in the presence of f (1) symmetry 
breaking terms has been discussed previously in Ref. [2T] . 
Next let us study the more interesting case where the 
dispersion e k has two degenerate minima at finite wave- 
vectors ±q. At the first sight it seems that in this case 
one can find solutions of the Gross-Pitaevskii equation 
(11 1 where only the modes with k = ±q condense, 



0fc = V^[<Wl +**,-»<!]■ 



(17) 



Keeping in mind that Vi = { a q)/vN and V-i = 
(at)/y/~N, we see that (V>f)* = V'-i- In real space the 
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condensate wave-function (171 corresponds to 



Setting ipf — ifj a _ 1 — ip, the corresponding condensate 
density is 



p 1 (r) = |^(r)| 2 =4|^| 2 cos 2 (q-r). 



(19) 



The important point is now that a condensate wave- 
function of this type does not solve the Gross-Pitaevskii 
equation (111, because the interaction terms couple the 



Fourier components with k = ±q to all other Fourier 
components involving arbitrary integer multiples nq of the 
fundamental wave-vector q, where n — 0, ±1, ±2, .... To 
see this more clearly, let us substitute the general ansatz 



(20) 



N Y, Sk tqn W 



n— — oo 



into Eq. (Ill where q n — nq. Setting the external wave- 
vector k = —q n in Eq. (Ill and defining r n = e_ g „ — /i 



and 7„ = r f- q „ (assuming again that 7^ is real), we obtain 
the discrete Gross-Pitaevskii equation, 



- y 



CT2 

n 2 



n±n 2 1J102 



u n,n 1 +n 2 +n 3 <- J nnin2ri3 <^ ni <^ n2 V 7 ^ J \^^) 



n 1 n 2 n 3 a 1 a 2 (T 3 



where 



T^<Tcri<T2 
' nn\n 2 



A(QtjCT, qn 1 o~i,q n2 o~2), 



(22) 



KnlnlZ = r A (q n a, q ni ai,g n2 a 2 , q n ^?,)- (23) 

The crucial point is n ow that if we assume on the right- 
hand side of Eq. (I2T]) that only the coefficients ip*'. with 
m = ±1 are finite, then we find after carrying out the sum 
that on the left-hand side all field components ip° with 
n = 0, ±1, ±2, ±3 must also be finite, so that the assump- 
tion that only the modes with wave- vector ±q condense is 
not self-consistent. For general interactions where all in- 
teraction coefficients V^nlnl an d U^nlnln\ are finite, the 
Fourier transform of a self-consistent solution of the Gross- 
Pitaevskii equation must therefore have finite weight for 
all integer multiples of q. Depending on the behavior of the 
interaction coefficients the spatial behavior of the conden- 
sate wave- function can look rather differently. In Sec. [3] we 
shall show that for YIG the cubic interaction coefficients 
Vnnlnl actually vanish identically, and that the behav- 
ior of the quartic coefficients Unn[n%\ IS sucn tn at the 
component ip±i of the condensate wave-function is much 
larger than the other components. In this case the spa- 
tial distribution of the condensate density is to a good 
approximation given by Eq. (19). On the other hand for 



some other types of interactions many Fourier components 
of the solution of the discrete Gross-Pitaevskii equation 
can have the same order of magnitude. In this case the 



5 - 



4|-0i | 2 



1 - 








0.5 



xq/2ir 



1 



Fig. 2. (Color online) This plot illustrates the fact that 
the condensate wave-function is more localized in real space 
if many Fourier components contribute. The dash-doted and 
the solid lines represent condensate densities associated with 
a condensate wave-function where the 5 odd Fourier modes 
Qi)Q3,Qs,q7, <79 are macroscopically occupied. For the dashed- 
dotted line we have assumed that all 5 modes have the same 
weight, while the dashed line corresponds to YIG with pump- 
ing parameter 71/ri = 3, see Table [l] For comparison, the 
dashed line is the condensate density pi (r) where only the 
modes with wave vectors ±q are macroscopically occupied, see 
Eq.FT9l 



condensate density is strongly localized at the sites of a 
one-dimensional lattice with spacing 27r/|q|. For example, 
if we assume that the first to Fourier components are fi- 
nite and equal, ?/>^ = ip/i/m for 1 < \n\ < to, then the 
condensate density is given by 



Pm(r) = 



2ip 



Y cos(q n ■ r) 



71=1 



(24) 



To obtain a properly normalized density, it is necessary to 
scale the order parameter as t/j/y/m. In Fig. we compare 
the single-component density p\(r) given in Eq. (19) with 



the corresponding multi-component density where the 5 
odd Fourier modes 91,93, 95, 57, <79 are macroscopically 
occupied. Obviously, in this case one can already observe 
a strong localization of the condensate at the sites of a 
one-dimensional lattice with spacing n/\q\. It is then ap- 
propriate to think of BEC as a condensation phenomenon 
in real space. In fact, the formation of the Bose condensate 
resembles in this case the phase transition from a liquid to 
a crystalline solid [9 . However, in the three-dimensional 
crystal formation problem the situation is more compli- 
cated because the Gaussian term in a Ginzburg-Landau 
theory exhibits a minimum on a surface in momentum 
space, and for the crystal structure the cubic term in the 
expansion of the Landau functional in powers of the den- 
sity play also an important role [TOlITT]. 



3 BEC of magnons in YIG 

It is well known [71IT7][2^1[2^1[271[2^pgil3l}l[3T] that the ef- 
fective magnon Hamiltonian for YIG can be cast into the 
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general form given in Eqs. (jlfHl- In the appendix we sum- 
marize the main steps and approximations in the deriva- 
tion of the magnon Hamiltonian for YIG from a real- 
istic spin Hamiltonian and give explicit expressions for 
the interaction vertices. If the samples have the shape 
of thin stripes and if an external magnetic field is ori- 
ented along the direction of the stripes (which we call 
the ^-direction), the energy dispersion e^ of the lowest 
magnon band indeed has two degenerate minima wave- 
vectors ±q = ±ge z . From the explicit expressions for the 
three-point v ertices for YIG given in Eqs. (A16 A17) and 
( A28a| A28d I we see that for this direction of q all thrcc- 
magnon interaction vertices V°nlnl defined in Eq. (22 1 
vanish identically, so that for the discussion of BEC in 
YIG one can omit the first term on the right-hand side of 
the discrete Gross-Pitaevskii equation (21 1. We can then 



construct self-consistent solutions of this equation involv- 
ing only Fourier components ip± n with n = (2j + l)ni, 
where j — 0, 1, 2, . . ., and n\ is the label of the lowest finite 
Fourier component. Because in experiments the samples 
are kept at room temperature, leading to a finite ther- 
mal magnon density at k = ±q, and energy transfer by 
pumping is mostly done to modes whose energy is less 
than e 2 qi we expect that the Fourier components i/'ii are 
dominant. In the following we therefore set m = 1 and 
consider solutions of the type 



E 

n odd 



S k.,q n 1pn 



(25) 



The infinite set of Fourier components ip±i, ip±3, ip±5, ■ ■ ■ 
is determined by setting n = ±1, ±3, ±5, ... in the discrete 
Gross-Pitaevskii equation (21), keeping in mind that for 

and 
d23b, 



BEC of magnons in YIG we should set V™*£ 
use the four-point vertices UZZ^Z 2 Z 3 defined via Eqs. 

— ^ — | |^^_^^_| I ti«l /i>2 '•'3 

(A18-A20), and (A29a-A29e). We have solved these equa- 
tions numerically by truncating the expansion ( |25| at 
some finite order m > \n\. For positive n non-trivial so- 
lutions can be obtained for 



7i > r i 



f 



(26) 



If ri is negative, we find solutions for arbitrary 71, includ- 
ing 71 = 0. As discussed in the appendix, to describe the 
stationary non-equilibrium state of the magnon gas in YIG 
under the influence of an external microwave field oscillat- 
ing with frequency ujq, one should re-define e q — > e q —ujo/2 
and use an appropriate chemical potential [i. It turns out 
that the Fourier coefficients ip° decay rapidly for large n, 
so that in practice it is not necessary to choose the cut- 
off to larger than 10 to obtain converged results. Typical 
numerical results are summarized in Table I and are rep- 
resented graphically in Fig. |3j 

Obviously, for YIG the first Fourier components tp^ 
of the condensate wave function are dominant, so that 
the spatial structure of the condensate density is to a 
good approximation given by p\ (r) given in Eq. ( 19 ) . This 



is consistent with the experiments by Demokritov and 
co-workers [TJUHMS] > wn0 observed a strong enhance- 
ment of the magnon distribution only for wave-vectors 



Table 1. Numerical results for |t/>i| = \ip\\ and the ratios 
l^nl / |^i| for different values of the dimensionless pumping 
parameters 71/ri for n > 0. Note that for r\ > there is 
no condensate if 71 /r\ < 1. The numbers have been obtained 
from the numerical solution of the discrete Gross-Pitaevskii 



equation (211 using realistic interaction parameters for YIG. 



By symmetry, for negative momenta we obtain identical re- 
sults. 

71/n 1^1 1 IWV'il IV'sMI IV'tMI IV^MI 



1.1 


0.530 


0.019 


0.002 








1.2 


0.757 


0.033 


0.005 


0.001 





1.5 


1.209 


0.055 


0.013 


0.003 


0.001 


3.0 


2.439 
\i>n\ 


0.085 


0.031 


0.012 


0.005 
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Fig. 3. (Color online) Absolute values \ip n \ = IV'nl °f the 
Fourier components of the order parameter for BEC in YIG for 
different values of the dimensionless ratio 71/7-1 = 7q/(e q — [l)- 
The order parameter has been obtained from the numerical 



solution of the discrete Gross-Pitaevskii equation (21 1, using 



the fact that for YIG the three-point vertices Vnn lnl. vanish. In 
this case the physically relevant solution of Eq. (21 1 has only 



odd Fourier components. Note the Fourier components t/>i are 
dominant and that the higher order Fourier components ip„ 
decay approximately exponentially as a function of n. 



±q. Note, however, that in principle also the higher or- 
der Fourier components are finite. In fact, with increasing 
amplitude of the pumping field the relative weight of the 
higher order Fourier components also increases. For ex- 
ample, from our numerical data shown in Table I we see 
that for 71 /n = 3 the amplitude of the third Fourier com- 
ponent ^±3 is approximately 8% of the amplitude of the 
dominant component ip±i- If the first Fourier components 
have the same order of magnitude and one uses realistic 
parameters for YIG, we estimate that the effective lattice 
spacing is 2n/\q\ ~ 2 • 10~ 7 m, which is roughly a factor 
of 100 larger than the spacing of the underlying Bravais 
lattice. 



4 Summary and conclusions 

In summary, we have considered the general problem of 
BEC in an interacting Bose gas whose energy dispersion 
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has two degenerate minima at finite wave- vectors ±q. Our 
main result is that for generic interactions the condensate 
wave-function has finite Fourier components ip n for all in- 
teger multiples q n = nq of the fundamental wave- vector q. 
For special interactions many Fourier components %l> n can 
have the same order of magnitude, so that in real space 
the condensate is strongly localized at the sites of a one- 
dimensional lattice. In this case there is a formal analogy 
between BEC and the liquid-solid transition. 

We have also used our theory to study the condensate 
wave-function of the condensed magnons in the magnetic 
insulator yttrium- iron garnet. In this case it is appropriate 
to think of magnon BEC as a condensation process in mo- 
mentum space, because the condensate wave-function is 
dominated by its leading Fourier components at ±q. How- 
ever, if the amplitude of the oscillating external microwave 
field is increased, then higher order Fourier components of 
the condensate wave-function can be populated. 

If one succeeds to extend the experiments which ob- 
serve directly the magnon densities to wavevectors k > q 
it should be possible to detect also these higher order com- 
ponents. Since the order of magnitude of the higher order 
Fourier components drops exponentially and derealiza- 
tion problems occur due to a finite slope of the dispersion 
away from the minimum, it is a rather challenging task to 
verify our prediction experimentally. 

Finally, it should be mentioned that quite recently 
Malomed et al. [33 studied the dynamics of BEC of 
magnons in YIG within an approach where coupled time- 
dependent Gross-Pitaevskii type of equations for the two 
components of the condensate wave-function correspond- 
ing to condensation at ±q are written down phenomeno- 
logically. Note, however, that this approach neglects pro- 
cesses which couple the dominant Fourier components t^±\ 
to the higher order Fourier components of the condensate 
wave- function. 
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APPENDIX: MAGNON-MAGNON INTER- 
ACTIONS IN YIG 

In this appendix we outline the derivation of an effective 
boson Hamiltonian of the form given in Eqs. (fl}[4| de- 
scribing the lowest magnon band of YIG, starting from 
the following time-dependent spin Hamiltonian, [7l|T6 ,17, 
|2"2"P3"P4] 



ij a0 



where a,0 = X,y,z label the three components of the spin 
operators Sf, and i,j = 1, ...,N enumerate the N sites 
t*j of a cubic lattice with lattice spacing a f=a 12.376 A. 
The exchange couplings Jy — J(ri — Tj) have the value 
J « 1.29K if Ti — Tj connect nearest neighbor sites and 
vanish otherwise. The Zeeman energy associated with a 
static external magnetic field H e is denoted by ho = \xH e , 
where /i = <?/is with the Bohr magneton given by \ib ■ Set- 
ting g — 2 we should work with an effective spin S « 14.2 
as d iscus sed in Refs. [71115]. The time-dependent part of 
Eq. ( Al ) represents the Zeeman energy induced by an ex- 
ternal microwave field oscillating with frequency cuq. The 
energy scale hi associated with the oscillating component 
of the magnetic field is assumed to be small compared 
with ho, so that both the static and the oscillating mag- 
netic field point into the direction of the macroscopic mag- 
netization, which we call the z-direction (parallel pump- 
ing). Finally, the matrix elements of the dipolar tensor 
D"f =D a P(r l -r J ) are 



Dff = (l-S ij ), 






STQ/9 



(A2) 



where ry = r.; 



Vj and Vij — fij/ 



Because the experimentally relevant YIG stripes are 
several thousand lattice spacings thick, we may assume 
that for magnetic fields oriented along the direction of the 
stripes the classical ground state is a saturated ferromag- 
net with all spins pointing in the direction of the external 
magnetic field. The components of the spin operators can 
then be expressed in terms of canonical boson operators 
hi and b] as follows 1231. 



S+ = V2S\ I 



sr = 




i*=^ 



'2S 



b\bjbi 
45 



b\b\bi 
AS 



O., 



S-b\bi, 



(A3a) 



(A3b) 

(A3c) 



where S+ = Sf+iSf and S t — Sf — iSf . Retaining terms 
up to fourth order in the boson operators we obtain 

H Y1G (t) = H Q (t) + H 2 (t) + H 3 + H± + OiS- 1 ' 2 ), (A4) 

where the boson-independent term and the term quadratic 
in the bosons are [71IT7] 



ffo(t) = — =-£[ J « + D%] 



-NSiho + hxCOS^ot)], 



H2(t)=J2 



Aijblbj 



B, 



\bibj 



b\b\ 



- [ho + hi cos(w t)] £ SI 



(Al) 



-/i 1 cos(w t)^&J& i , 



(A5) 



(A6) 
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with coefficients given by 



The interaction parts can be written as 

#3 = —7= E <5fci+fc 2 +fc 3 ,0 



N 



B, 



-S 



S 



5- Vr- ij ij 



-\n xx - 2in xy - n vv ] 



(A7) 



(A8) 



The cubic contribution H3 to the boson Hamiltonian is of 
order yS and involves only the dipolar tensor, 



H A = 



1 



1 

N 



I [rl%bt x b 2 b 3 + r^bljUh 

E '*- 



2! 



(A14) 



fcl.../C4 



-fc 2 + fc3+fc4,0 
1 



1 



^bbbb it . t 



(2!) 2 



rs&bijizhh 



~r^ 4 bUb 2 b 3 h + ^r^bUbUbi^ 



H a 



l/lEtt 



v 



■X 



Dff+iD^) lb]b]b< + ~b]b]b 



Ut. 



1 



.ti.ti 



■ j J T 4 



h.c. 



(A9) 



(A15) 
where the properly symmetrized three-point vertices are 



r bbb 
1 1:23 



^ " ^ + (fe -> fe) 



The quartic part of the boson Hamiltonian can be written 

as 



(D?-iD?) 



H i = ~2 E Ji J [ ntU J ~ 2 \fA h J b 3 + h 

+ ^ E ( D « + c « ) h^ + 1 ( 6 ! 6 M + h - c -) 



(A16) 

^sf a = rei)*, (ait) 

and the symmetrized four-point vertices are 



r bbbb 
1 12:34 



{D%-2iD%-D%)b\b t b i b j +h.c. 



Jkx+k.3 + Jk 2 +k 3 + ^fei+fc 4 + Jk 2 + k A 

r> zz 4- n zz 4- n zz 4- n zz 



(A10) 



-£(J fci -2D|*) 



where we have abbreviated rij = 6^ 6^ . Next, we Fourier 
transform the Hamiltonian to momentum space, setting 



\Y 



h = 



1 V^ ikr ih 



(All) 



r bbbb _ ■*• I nil • n a;j/ n i/j/ 
-* 1;234 - 4 |_-^fc 2 _ / * i7 fc 2 ~ U k 2 

+ (fc 2 ~» ^3) + (&2 ""► fc 4) 

r bfc6b _ / r bbbb \ * 
1 123;4 — 1/ 4;123j ' 



(A18) 

(A19) 
(A20) 



where for simplicity we impose periodic boundary condi- 
tions in all directions. A more accurate calculation should 
take into account the finite extend in the direction where 
the experimentally relevant samples have the smallest ex- 
tension (which we call the x-direction) [7 . For our purpose 
it is sufficient to impose periodic boundary conditions in 
all directions, which amounts to approximating the eigen- 
functions of the exchange matrix J^- by plane waves. The 
lowest magnon band is then obtained by simply setting 
k x — 0. In Ref. [7] we have shown that this uniform mode 
approximation reproduces the qualitative features of the 
dispersion of the lowest magnon mode rather well. In mo- 
mentum space, the quadratic part Hi{b) of our bosonized 
Hamiltonian becomes 



The Fourier transforms of the exchange and dipolar cou- 
plings are defined by 



ik-rij j 



'* = £ 

i 



(A21) 
(A22) 



Finally, we use a Bogoliubov transformation to diagonalize 
the time- independent part of Hi{t), 



Hi 



E 



A k V k bk 



Hk 



(bibt h 



b~kbk 



bk 



-k 



Uk ~Vk \ I a k 

-v* k u k ) \a}_ k 



where 



+hi cos(u t) ^2 b k bk, 



(A12) 



_ I A k + e k B k A k - e k 



with 



A k = E e "* r *^j> B * = E e-' fery %- (A13) 



and 



e k = JAl - \B h \*. 



(A23) 



(A24) 



(A25) 
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After this transformation the quadratic part of the Hamil- 
tonian reads [Tfl[25] 



ff a («) = £ 



e k a k a k 

h 1 cos(cl> £) 
cos(woi) 



Cfe 



A k 



E 

k 



A k t 
— a k a k 



k 



-k 



where 



Ik = 



pal 



Ak — £fc 
2e k 

7fc a -fcOfc 



2e fc 



(A26) 



(A27) 



Substituting the Bogoliubov transformation JA23 ) into 



the expressions for H 3 and H 4 given in Eqs. (A14 A15), 
we arrive at expressions given in Eqs. ([3] [41, with the cubic 
vertices explicitly given by 



nbbb 



bbb 



^bbb 



+r™? 3 v lV2 u 3 + ri^v 2 v 3Ul + r^%v lV3 u 2 , 



(A28a) 



A a ;2?3 = JS«i«2«3 + rl%v lV2 u 3 + r^ 2Vl v 3 u 2 



-r^ l v a v 2 v 1 - r^ 3 v 2Ul u 3 - I^Witia, 



r-iaaa { j^aaa\* 

1 12;3 ~~ I 1 3;2lJ > 

7-iaaa ( -paaa\* 

1 123 — 1/ 123 J ' 



(A28c) 
(A28d) 



The vertices appearing in the quartic part H 4 of the 
Hamiltonian in the Bogoliubov basis are (see Eq. Q) 

/™ = r|^tti«2W3t;4 + r^UlU 3 V2VA 

+/f 4 %UlU 4 ^3 + ^iWs^ + r^%U 2 U 4 V lV3 

+rgf 2 U 3 U4«l«2 " 4 6 l b 23WlU2U 3 V4 - rg&uniatuva 

-rg^UxUatUWS - r^«2« 3 «4Vl - rjgiui«2t>3V4 

-r|i? 2 U2Wl«3«4 - ^4 b 3 W 3«l«2«4 " A^^l^S, 

(A29a) 
A a ;234 = -I%&U2ViVaV4 - rjg&ujvufct* 

-rl\%uiviv 2 v 3 - r 3 ^ 2 u 3 u 4: uiv 2 - r 2 ^ 3 u 2 u 4 uiv 3 

-r*l%u 2 u 3Ul v 4 + fg&ttitfcttau* + r]^iu3«2«i«4 

+rUuiu4«2«i«3 + rl%u 4 u 3 v 1 v 2 + r^ A u lUl v 2 v 3 
+rHl%U 3Ul V2V4 + r^u 2 u lV3 v 4 + rjg^t^vi, 

(A29b) 



bbbb 



A°2°34 = A2;34 U 1 U 2«3W4 + A3 ; 42«l«*«3«a 



+r 1 fe 4 & . f 3 h 2 u lU3 w4«2 + r^V^^i + r| 4 fc f 3 V 2 u 3 « 4 f 1 
+r^% Vl v 2 v 3 v 4 - rl%u 3 v 2VlVi - rl%u 4 v 2 v lV3 

-r|£&U2U3U4Vl " ^342^3^2 - r^i«2«3«4Vl 

-r , 1 6 ^ 6 2 Ul-y 3 U4U2 - A 6 24 6 3 u l'"2W3«4 ~ A 6 23 b 4 M l u 2M4«3, 



134;2 l 



(A29c) 



and 



^aaaa 
1234 



r^aaaa 
1 1234 1 



j^aaaa / j^aaaa 

1 4;321 — \ l 123;4 



(A29d) 
(A29e) 



For nearest neighbor coupling on a cubic lattice with spac- 
ing a the Fourier transform of the exchange coupling ap- 
pearing in the above expressions is 

J k = 2 J[cos(k x a) + cos(k y a) + cos(fc z a)] . (A30) 

The Fourier transform of the dipolar tensor is more com- 
plicated. For a thin YIG film the minimum of the disper- 
sion is at ±q = (0,0, ±q), so that in this work we only 



need D? as a function of k? for fc T 



ty =0. For a film 
with thickness d 3> a we then obtain in uniform mode 
approximation [7], 



(A31a) 

(A31b) 
(A31c) 



U k z 


47T/X 2 

a 3 


l~ h '_ 


f 


di: 


47T/i 2 

3a 3 ' 




di: 


Airy 2 
a 3 


~3 +fk ° 


K 


= di: = 


D Z = ° 


1 



(A28b) where the form factor f kz is given by |S1[7] 

1 _ e -[fe»M 



fk„ — 



\k z \d 



(A31d) 



(A32) 



For the purpose of studying the phenomenon of para- 
metric resonance [2l[2^1[27lfe8l[2"9l l3l? l [3T ] . one usually sim- 
plifies the Hamiltonian ( A26) by dropping the second line. 



involving the combination cos(ujQt)A k a k a k ; moreover in 
the last line one substitutes 



7^ cos(w t) 



^e-^ , 7 *cosM)^^ 0t . 

(A33) 
Although the validity of this approximation in the con- 
text of YIG is questionable [IT] (see also Ref. [32]), let us 
assume here that it correcly describes at least some as- 
pects of the experiments 1 , 2 , 3, 4 , 5j. Our quadratic boson 
Hamiltonian is then approximated by 



H 2 (t) 



^ e k a k a k 

k 



2^h e 



- iuJot a{a\ 



* ILJnt 

7fc e a-kO-k 



(A34) 



where we have dropped the constant terms. The explicit 
time-dependence may now be removed by a canonical 
transformation to the rotating reference frame, 



ak 



e^° t/2 a k , 



st 



e-^ ot/2 a k , 



(A35) 



so that the transformed quadratic part of our Hamiltonian 

is 



Ho 



E 

k 



h 


a k a k 


+ 


7fe 

2 


4 a - 


k 


+ 


# 

7fc~ 

i a - 


-fcOfc 



(A36) 
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where ik = e& —ojq/2. If we re-define again a^ — \ a,k, If. 
efc, we arrive at the quadratic Hamiltonian in Eq. (pi). 
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